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ABSTRACT 

An ion beam can destabilize Alfven/ion-cyclotron waves and magnetosonic/whistler waves if the beam speed 
is sufficiently large. Numerical solutions of the hot-plasma dispersion relation have previously shown that 
the minimum beam speed required to excite such instabilities is significantly smaller for oblique modes with 
k X Bo y than for parallel-propagating modes with k x Bq = 0, where k is the wavevector and Bo is the 
background magnetic field. In this paper, we explain this difference within the framework of quasilinear theory, 
focusing on low-/? plasmas. We begin by deriving, in the cold-plasma approximation, the dispersion relation 
and polarization properties of both oblique and parallel-propagating waves in the presence of an ion beam. We 
then show how the instability thresholds of the different wave branches can be deduced from the wave-particle 
resonance condition, the conservation of particle energy in the wave frame, the sign (positive or negative) of 
the wave energy, and the wave polarization. We also provide a graphical description of the different conditions 
under which Landau resonance and cyclotron resonance destabilize Alfven/ion-cyclotron waves in the presence 
of an ion beam. We draw upon our results to discuss the types of instabilities that may limit the differential 
flow of alpha particles in the solar wind. 

Keywords: magnetohydrodynamics (MHD) -solar wind - Sun: corona - turbulence - waves 



L INTRODUCTION 

The solar wind is a dilute plasma consisting of protons, 
electrons, and other ionic species. Among these other ions, 
the alpha particles play a special role due to their relatively 
high abundance. In the fast solar wind, for example, alpha 
particles compri se ^ 1 5-20% of the solar-wind mass den- 
sity dBame et alJll977h . which makes them dynamically and 
thermodynamically important. 

In situ measurements have shown that the alpha parti- 
cle velocity distribution is far from collisionally relaxed. 
For example, alpha particles generally have higher tem- 
peratures and outflow velocities than the protons. Also, 
like protons, alpha particles exhibit temperature anisotropy 
with Tj_ y 7]|, where T± (7||) measures the speed of ther- 
mal motions perpendicular (parallel) to the local mag- 
netic field (Marsch et al. 1982; von Steiger et al. 19951 



Neugebauer et al.iil996c iReisenfeld et al.ii2001t iMaruca et al.l 



When measured in the proton frame, the beam speeds 
of alpha particles Ui seen in situ at heliocentric distances r ex- 
ceeding 0.3 AU is in most cases limited to a value of order 
the local Alfven speed va- At these heliocentric distances, the 
values of both va and Ui decrease with increasing r while the 
proton speed is nearly constant. This means that an ongoing 
deceleration mechanism has to act on the alpha particles and 

regulate their drif t speeds in in terplaneta ry space. 

Recent studies jKasper et al. 2008; B ourouaine et alj|201 Ih 
have shown that Ui is anti-correlated with the collisional age 
rVco\/U of the solar wind, where i^coi is the Coulomb colli- 
sion frequency and U is the proton outflow velocity. Un- 
der cond itions of rare colUsions i n the measured solar wind 
interval, iBourouaine et al.l (1201 Ih found an upper threshold 
of about 1.2va f or Uj in Helios data at r 'y 0.7 AU. Sev- 
eral au t hors (e.g..llsenberg & Hollweg|1983l;lGomberoff et al.l 
ll996btlKasperet al.ll2008i;lAranedaetal.ll2009h have argued 
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that wave-particle intera ctions would lead to an upper limit 
of f/i ^ Va. The study bv ' Bourouaine et aTl ( 1201 Ih also found 
that solar-wind streams with Ui near this upper limit show ex- 
cess magnetic fluctuation power at wavenumbers of up to the 
inverse local proton inertial length. 

In situ observations suggest that the solar wind hosts a 
number of different plasma wave modes, including Alfven 
waves (Belcher & Davis 1971; Tu & Mai'sch 1995), ion- 
cyclotron waves (Jianetal. 2 010t iHe et al.l 1201 2D , whistler 
waves (Podesta & Garv. .201 li) , and kinetic Alfven waves 
dBale et alJ 120051: iSalem et alj 120121) . The properties of 
these waves depend on the properties of the plasma in 
which they propagate. One of these background prop- 
erties is the composition of the plasma. Adding other 
ion species or relative drifts between the species mod- 
ifies the w a ves s ig nificantly ( Gomberoff & Elgueta 199ll 
Gnavi e t al.'' '19961; 'Gomberoff et al. 1996a; PeiTone et aD 



2011 



2011 



Marsch & Verscharen 201 1; Verscharen & Marsch] 
The above-mentioned alpha-particle beam is a typ- 
ical example of such a drift. Apart from the modifica- 
tion of dispersion properties, the beam provides a source of 
free energy that can lead to the growth of electromagnetic 
instabilities, as shown in analytical and numerical studies 
(e.g., Montgomerv et al."1976VGarvl ll993l;lGarv et al Il2000bl 
Li&Habbal 2000; Lu et al. 20061 120091: iLi & Lul l20ioh . 
These drift instabilities decelerate the drifting particle species 
to a stable value, thereby consuming the free energy source. 
These treatments discuss different types of instabilities and 
determine the threshold value of the relative drift at which 
the plasma becomes unstable. Some of the instabilities dis- 
cussed include a circularly polarized parallel magnetosonic 
mode, an oblique magnetosonic mode, and two Alfvenic 
modes. The magnetosonic mode reaches its maximum growth 
rate at parallel propagation, at which th e (previously stud- 
ied) Alfvenic instabilities are suppressed ([Montgomery et al] 
119761) . iGarv et al.l (I2000bl) have shown that the threshold drift 
speed for oblique Alfvenic instabilities is significantly smaller 
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than for the parallel magnetosonic instability — Ui'>j 0.8va 
instead of Ui > 1.7va, at least for situations in which the elec- 
tron temperature is fairly large (four times the proton temper- 
ature in their study). 

The intention of the present article is to give an analytical 
explanation for this behavior within the framework of quasi - 
linear theory. Our primary objective is to elucidate the phy sics 
of these instabilities and not to refine or revise previous re- 
sults on the values of these instability thresholds. We base 
our analysis on the cold-plasma dispersion relation, which we 
expect provides a reasonably accurate description of the real 
parts of the wave frequencies when /3 <C 1, where /3 = ^TTp/B^ 
is the ratio of the plasma pressure p to the magnetic pressure. 
However, our use of the cold-plasma dispersion relation limits 
our analysis to the solar corona and the subset of solar-wind 
streams farther from the Sun that satisfy /3 <C 1 . Since most 
of the solar wind near 1 AU satisfy /? ^ 1, our results do not 
describe the bulk of near-Earth solar wind. The remainder of 
this paper is organized as follows. In Section |2] we derive 
the cold-plasma dispersion relation for obliquely propagating 
plasma waves in the presence of an ion beam. We use the so- 
lutions of this dispersion relation to evaluate the wave polar- 
ization and wave energy. In Sections[3]and|4]we describe how 
the instability thresholds can be determined and understood in 
terms of (1) the resonance conditions for wave-particle inter- 
actions in quasilinear theory, (2) the conservation of resonant- 
particle kinetic energy in the reference frame moving with the 
wave along the background magnetic field, (3) the wave po- 
larization, and (4) the sign (positive or negative) of the wave 
energy. In Section |5] we summarize our findings and discuss 
their application to alpha-particle beams in the solar wind. 

We note that we do not address the details of the mech- 
anism(s) that generate ion beams in the solar wind, which 
could be cyclotron-res ona nt wav e- particle interactions 
(e.g., McKenzie & MarscH Il982t iHollweg & Isenberg 
|2002; Ofman et al. 2002 ), inhomogeneity effects with low- 
frequ ency waves (McKenzie et alj[l979l : llsenberg & HoUwegl 

119821). stochasti c he ating by low-frequency waves 

jMcChesnev et al."198r; 'Johnson & Cheng"2001'; 'Chen et al 



2001; Chastonetal. 2004; Chandran 2010; Chandran et 
20101) . and/or some alternative process(es). Also, we con 



all 



centrate on instabilities driven by resonant particles and not 
on the beam instabilities that a re seen in solutions to the 
cold-plasma dispersion relation (Gomb eroff et al.|[l996al) . in 
part because these "cold" instabilities have been found to be 
stabilized by the presence of other waves (lGomberoffll2"00l 
lAraneda & Gomberoff„200l . 

2. WAVE PROPERTIES IN A COLD PLASMA 

We consider a uniform, cold plasma consisting of three par- 
ticle species: protons, electrons, and a population of beam 
ions. We assume that there is a uniform background mag- 
netic field Bo = Boe,, and that the electrons and beam ions 
flow with respect to the protons at average velocities f/eCj and 
f/iC,, respectively, with Ui> and Bq > OQ We carry out all 
our calculations in the proton rest frame. We assume that on 
average the plasma is charge neutral, with 



npg + Mi^i- He^ = 0, 



(1) 



and carries (on average) zero current, 

Mi^it/i-Meet/e =0, 



(2) 



where Wp, tii, and rig are the proton, beam-ion, and electron 
number densities, and e and q-, are the proton and beam-ion 
charges. 

2.1. Dispersion Relation 

We write the Fourier transform of the electric field in terms 
of its Cartesian components E^x, E^y, and /ij^ EI By linearizing 
Maxwell's equations and the momentum equations for each 
particle species, we obtain an equation of the form 



kc /kc 

X I — xEk 



+ £E* = X>Et=0, 



(3) 



where k = k±ex + k^^e. is the wavevector, uj is the frequency, 




(4) 



is the dielecti-ic tensor, S = (R + L)/2, D = (R-L)/2, and the 
quantities 



R=l 



■E 



\ UJ 



and 



uj — k^^Us 
uj — k^\Us + ^ls 

LO — k\\Us 

u}-ki\Us-fis 



(5) 



(6) 



extend the notation of 'Stij^l d 19921) to account for the rela- 
tive drift between particle species. The plasma frequency of 
species s is given by lu^^ = Airnsq^Jnis. The indices _L and | 
refer to the directions perpendicular and parallel to the back- 
ground magnetic field. The gyrofrequency of the correspond- 
ing particle species is denoted fij = qsB/(msc) and is negative 
for electrons. The terms 



Us 



uj — k\\Us 



and 



M-- 



(c.-^l|f/,)'-n? 



(7) 



(8) 



represent new components of the dielectric tensor that are ab- 
sent when the beam is absent, and 



P=l 



■E 



UJ 



{uj-k\\U,) {uj-k\\U,) -n] 



(9) 



The term "1" in the expressions for R, L, and P stems from 
the displacement current and can be neglected in our study. 
To obtain the dispersion relation, we set 



det2:' = 0. 



(10) 



We solve Equation ( fTOl i in two steps. First, we obtain ap- 
proximate solutions to the dispersion relation in the limit that 
me 0, where nig is the electron mass. Second, we use this 



^ In the anti-parallel case (U; ■ Bq < 0), the same theoretical arguments 
apply, and the same instabilities are found. Some labels, however, have to be 
changed in this case. 



^ The quantity is t he spatial Fourier transfomi of the electiic field, as 
defined in Equation {24). To obtain the dispersion relation for linear waves, 
we take Ei. cx: e^'^' . 
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approximate solution as the starting point for a Newton's- 
method solution to the full cold-plasma dispersion relation, 
taking into account the finite value of m^/m^. 
In the limit 0, Equations (|5]l through dH) become 



rip (r^p+cjj 



rip (r!p- 



and 



2rii(f7i+w-/t||t/i) 

2rii (Oi-w + ;t||t/i) 
uj-k\\JJi 



(w-fc||f/i) 

(w-A;||t/i) 



(11) 

(12) 
(13) 

(14) 



As nif. — > 0, P — > oo, and only those terms in detl? that contain 
a factor of P need to be retained when solving Equation (fTOl l. 
After factoring out this factor of P from Equation (fTOl l. we ob- 
tain the (m^, — > 0) dispersion relation for Alfven/ion-cyclotron 
and fast/whistler waves. 



R + L 



[l+cos^0)+RL = Q, (15) 



where 6 is the angle between k and Bq. When 9 = 0, Equa- 
tion dTsl l becomes 



= 0. 



(16) 



Each of the equations (kc/uj)^-L = and {kc/ijj)^-R = leads 
to a third-order polynomial equation for lo. When the beam 
ions are alpha particles, the equation (kc/uj)^—L = i s equiv- 
alent to Equation (2) of iGomberoff & El^uetal (Il991h . When 
9 ^Q, Equation (fTST i leads to a sixth-order polynomial equa- 
tion for Lu, which w e solve numerically using the Laguerre 
method dPress et al.lll992i) . 

As mentioned above, we use our solution of Equation ( fTSl l 
as the starting point for a Newton's-method solution of Equa- 
tion ( [Tol l taking into account the finite electron mass. This 
second step to our solution of Equation ( fTOl l leads to only mi- 
nor refinements to the wave frequency. However, it enables 
us to calculate polarization variables for each mode such as 
Ekz/Ekx, which are relevant to our discussion of Landau damp- 
ing in the sections to follow. 

For the remainder of this paper, we use the full solution of 
the cold-plasma dispersion relation with finite nif. rather than 
working in the zero-mg limit. We also from here on restrict 
ourselves to the case in which the beam ions are alpha parti- 
cles. We define 

m 

V=- (17) 
rip 



and 



Bo 



(18) 



In Figure [T] we plot Re w = for each of the six roots to 
Equation ( fTSl ) for three different choices of rj, U\, and 9. We 
focus on four of these branches, which we label according to 



their asymptotic behaviors at large k\\ . Branch A is the proton- 
cyclotron branch, for which lu ^ ftp as ^|| — > oo. Branch B 
is the fast-magnetosonic/whistler branch, which behaves like 
a fast magnetosonic wave at kv^/^p < 1 and a whistler wave 
at kv^/^p > 1. Branch C"^ is the forward-propagating alpha- 
cyclotron branch, for which uj k^^ U\ + Vli at large ^|| . In the 
frame of the alpha particles, the Doppler-shifted frequency 
of this wave Lo-k]fPi approaches ili as ^|| oo. Branch C~ 
is what we call the "counter-propagating" alpha-cyclotron 
branch. In the alpha-particle frame, this wave propagates in 
the -z direction and has a Doppler-shifted frequency Ljj-k\\Ui 
that approaches -fii as k\\ oo. We note that for each case 
shown, Imw = 7 = except at the intersections of different 
branches of the dispersion relation, where two branches have 
equal Wr- In the right panel of Figure [T] this occurs, e.g., at 
k\\VA/^p « 1.1 and cl>i- w O.Tfip, where 7 ~ ±0.02ilp. Both 
solutions have the same absolute value of 7 but with oppo- 
site signs. Other modes can also merge in certain wavenum- 
ber ranges for different parameter choices. The magnetosonic 
mode B and the counter-propagating alpha-cyclotron mode 
C~ merge at higher drift speeds and show also a non-zero 
growth rate in the wavenumber range where they have equal 
Wi . This unstable behavior of the cold-plasma dis persion rela- 
tion is known as the ion-ion r esonant instabihty (iGnavi et alJ 
[T996HGomberoff et al.l[T996ah . 

2.2. Wave Polarization 

The sense of circular or elliptical polarization of the waves 
can be defined in terms of the quantity 



V 



_ \EkA-\Ek,\\ ^ \A-i\-\A + i\ 
~ \Ek,\ + \Eu\ |^-/| + |^+/r 



(19) 



where Ekj = {Ekx-iEk,)/ \/2 and Eks = {Ekx+iEky)/y2. With 
the use of the x and y components of Equation we can 
write 

{S - n] - n^){J + rixu^) + iDM 



A: 



{S-nl)M-iD(J+nxn.) 



(20) 



where = {kc/Ld)sm9 and «, = (kc/uj)cos9. When 6' = 0, we 
can use just the x component of Equation ([3]) to obtain the 
simpler expression 

Following the standard convention, we say that a wave is left 
(right) circularly polarized if the electric field vector in the x-y 
plane rotates in the same sense as the cyclotron motion of a 
proton (electron). The connection between left/right circular 
polarization and the value of V depends upon the sign of the 
real part of the wave frequency Uy. When ujy > 0, left-hand 
circular polarization corresponds to V = -I, whereas right- 
hand circular polarization corresponds to T-" = 1. When Ui- < 0, 
it is the reverse: left-hand circular polarization corresponds 
to = 1 and right-hand circular polarization corresponds to 
V = -I. In general, waves are not circularly polarized, but 
elliptically polarized with P /±1. Linear polarization corre- 
sponds toV = 0. The value of V is shown in Figure|2]for the 
case of oblique propagation. 

It is instructive to consider the sense of polarization of the 
wave branches A, B, and in Figure [T] when 9 = 0. It 
can be seen from Equations (fTST i and fTH that V = ±1 when 
9 = 0, so that these wave branches are circularly polarized. 
More specifically, branch B is right-circularly polarized, and 
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Figure 1. Solutions to Equation ^3) for an alpha-particle beam. The letters are explained in the text. For the left panel, rj - 
middle panel, rj = 0.05, 9 = 1.5°, and U, = I.Ii'a. For the right panel, -q = 0.05, 9 = 45°, and U, = I.Ii'a. 

modes to have Ek: = 0. 
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Figure 2. The dispersion relation is shown for 9 = 45°, r) = 0.05, and Ui = 
1.2va. The color coding represents the polarization function V. All modes 
in this diagram are elliptically polarized at oblique propagation. 

branches A and C"^ are left-circularly polarized. Branch C~ is 
also left circularly polarized when Ui = 0, or, if Ui > 0, when 
ojr < 0. However, when Ui>0 and cjr > 0, branch C~ is right- 
circularly polarized. It is important to recall, however, that 
we have defined the sense of polarization in the rest frame of 
the protons. If we were to switch into the rest frame of the 
alpha particles, the Doppler-shifted wave frequency LUy - U 
of branch C~ would always be negative, and in this reference 
frame branch C~ would appear to be left circularly polarized. 

The longitudinal polarization is important for Landau- 
resonant interactions. We define the longitudinal polarization 
as 



Ekz 



(22) 



The right-hand side of Equation (122b can be evaluated with the 
use of the first two lines of Equation (O. It is found to be al- 
ways much less than one for the treated modes. Its maximum 
value in the treated frequency and wavenumber range is of the 
order of Vr < 10"^. In parallel propagation, the determination 
of the parallel component of the electric field is not possible 
anymore via the first two components of Equation (|3) since 
the dielectric tensor becomes block diagonal. The only paral- 
lel mode with Ek, ^ is the P = Q mode, which has very high 
frequencies for the plasma conditions that we consider, and 
which is excluded here. Therefore, we can assume all parallel 



2.3. Wave Energy 

The wave energy density per unit volume in k space, de- 
noted W, is given by 

d 



(2TryW: 



1 



duj 



(23) 



where (B^) is the Fourier transform of the fluctuating elec- 
tric (magnetic) field, and eh = (e+e^)/2is the hermitian part of 
the dielectric tensor (Stix 1992, see also the discussion follow- 
ing Equation 27 of lChandran et aO2()I0l) . The dagger symbol 
denotes the adjoint matrix. The factor of (In)^ on the left- 
hand side of Equation ( |23] ) does not appear in Equation (20) 
of Chapter 4 of lStixl d 19920 . This difference arises because we 
define 



Et 



E(x,f)e"'''"d^ji:, 



whereas [Stixl(ll992h defines 

Et = (27r)"^/2 / E(x,f)e"''' M^x. 



(24) 



(25) 



This energy density describes both the energy density of the 
electromagnetic wave fields and the kinetic energy density of 
the particles under the assumption of a slow change of the 
wave amplitude with time. 

It was found by Kadomt sev et alj n965h that the term con- 
taining the derivative of the dielectric tensor in the expres- 
sion for the wave energy can be negative and even lead to a 
situation in which waves have negative energy, meaning that 
W < 0. In this case, the averaged kinetic energy of the parti- 
cles plus the energy of the electromagnetic field is larger when 
the wave is absent tha n when it is present (Da wsonllI96II; 
iHollweg & V515 [T97Th . The behavior of negative-energy 
waves under the influence of dissipation is counter-intuitive. 
Dissipation, which removes energy from the fluctuations and 
converts it into heat, causes negative-energy wa ves to grow 
(iHaU & Heckrottel[T96a iHollweg & Voll3ll97Th . Similarly, 
if a negative-energy wave loses energy through a paramet- 
ric instability, its own amplitude grows, increasing the rate 
of this instability. This is the so-called explosive instability 
dOraev skv 1984). We find that the counter-propagating alpha- 
cyclotron wave C~ has negative energy for all lj^ > 0. This sit- 
uation is depicted in Figure |3] Negative-energy waves in the 
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Figure 3. Orange (black) color curves indicate positive (negative) values of 
the wave energy W for the different solutions to the dispersion relation for 
8 = 45°, »7 = 0.05, and l/j = I.Iva- The blue-dotted line indicates a;, =k^^Ui. 



presence of an ion beam have been investigated previously by 
iGnavi et al.l([T996h . 

In all cases that we have investigated so far, when the 
alpha-cyclotron branch merges with the proton-cyclotron or 
the magnetosonic mode, the resulting branches have negative 
energy throughout the finite interval of in which the two 
branches have the same value of Wj. This situation can be seen 
in Figure [3] where the forward propagating proton-cyclotron 
wave (branch A) and the alpha-cyclotron wave (branch C~) 
that propagates backwards in the frame of the alpha particles 
merge. 

3. QUASILINEAR THEORY 

Quasilinear theory provides an approximate, analytical 
description of wave-particle interactions in a collisionless 
plasma. One of the central assumptions of this theory is 
that the fluctuation amplitudes are sufficiently small that they 
can be described as a superposition of linear waves, and a 
particle's orbit can be approximated as an unperturbed helix 
aligned with the background magnetic field. In addition, the 
theory assumes that the growth or damping rates of the waves 
are much smaller than the linear wave frequencies. The start- 
ing point of this theory is the Vlasov equation. 



dfs 



— -I-V-— -I- — 
at ox 



E-I--VXB 



■0, (26) 



which describes the evolution of the distribution function 
/5(x, v,0 of particle species s with mass m, and charge un- 
der the influence of the electric field E and magnetic field B. 
In a constant homogeneous background magnetic field (Bq = 
BoCj), a perturbation analysis leads to the following relation 
for the slow temporal evolution of the distribution function: 



^= lim y 



dt v^co ^-^ Sn^m^ I Vv± 
«=— 00 ^ 



1 



Gv^(5(a;i:r — ^11 V|| —nns)x 
x\4,„^kfGf,d\ (27) 



where 



^IP'II 



d kiiv±_ d 
■ + ■ 



and 



= --j= [Ek,e''^J„+,{a,) + Eu,\e-'*J„.,{a,)\ 

+ ^£fe7„((T,) (29) 

jKennel & EngelmanrfI966l:lStixll992l) . The azimuthal angle 
of the wavevector k is defined as (p. The parallel and perpen- 
dicular components of the velocity in cylindrical coordinates 
are denoted V|| and v^, respectively. The real part of the fre- 
quency that is a solution to the dispersion relation is denoted 
LOiij. The amplitude function il),,^ contains the Bessel func- 
tion of nth order J„ with the argument CTi = k±y±/Vls. The 
quantity V is a volume that is related to our Fourier transform 
convention. We treat the plasma as infinite and homogeneous. 
In order to make the Fourier transform of a function /i(x) con- 
verge, we first multiply this function by a window function 
W that is 1 inside a cubical region of volume V centered on 
the origin and outside this volume. We then take the limit 

of y 00. 

Equation dZTl l shows that, in the framework of quasilin- 
ear theory, waves cause particles to diffuse in velocity space. 
Because of the delta function in Equation ( |27] |, particles of 
species s with velocity v are affected by waves with wavevec- 
tor k and frequency ujki only if the waves and particles satisfy 
the resonance condition 



(30) 



where n is any integer Equation (l30l l is called the Landau- 
resonance condition when n = and the cyclotron-resonance 
condition when n ^ 0. If there are only waves at a single k 
and Wij , then the diffusive flux of particles in velocity space is 
nonzero only within narrow bands of V|| that correspond to the 
solutions of Equation ( l30l l. Within those bands, the diffusive 
flux of particles is tangent to semicircles in the V||-v_l plane 
defined by the equation 



(v|| - Vph)^ + v^^ = constant, 



(31) 



(28) 



where Vph = w^r/^ii is the parallel phase velocity of the waves 
(Kennel & Engelmann 1966). When a spectrum of waves is 
present, velocity-space diffusion is not necessarily limited to 
narrow bands of V|| values, and particles can instead diffuse 
throughout an extended region in velocity space. For exam- 
ple, energetic protons undergoing cyclotron interactions with 
a broad spectrum of non-dispersive (A;|| <C ^p/va) Alfven 
waves with dispersion relation ojkr = ^||Va scatter along ex- 
tended semi-circular arcs centered on vaC^. In contrast, when 
wave-particle interactions are limited to the Landau reso- 
nance, particles can diffuse over an extended range in velocity 
space only if the spectrum of waves has a continuous range of 
Vph values, because the Landau resonance only arises when 
V|| = Vph. At each resonant value of V|| , Landau-resonant inter- 
actions lead to diffusion only in V[|, not in v^, and so the ex- 
tended diffusion paths for Landau-resonant interactions with a 
spectrum of waves are horizontal lines in the V||-vj^ plane. For 
our purposes, it is sufficient to consider the localized velocity- 
space diffusion that results from waves at a single k and ujh-. 
In this case. Equation dSTT l determines the vector direction of 
the diffusive flux, up to a 180 degree ambiguity. This ambi- 
guity can be resolved by noting that the net diffusive particle 
flux is directed from regions of higher particle concentration 
to regions of lower particle concentration in velocity space. 
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Figure 4. Dispersion relation (black, solid) and cyclotron-resonance condi- 
tion (red, dashed) for = 45°, f/j = 1.2va, and rj = 0.05. The blue dotted line 
represents ujy = iyt/i. The triple-dotted gray line corresponds to = 0. The 

spread of the resonance lines corresponds to a parallel ion beta of /9||i ~ 10"* 
for this choice of rj. 



Bessel functions have the property J„(Q) = S„x), which shows 
that, for "slab" waves with k± = 0, resonances occur only 
when n = -1, 0, or 1 . We assume that 



cri< 1, 



set 



(32) 



(33) 



and hence ignore resonances with \n\ > 1. Equation ( |32] | 
is satisfied for Alfven/ion-cyclotron waves and magne- 
tosonic/whistler waves provided < ftp, /3 ^ 1, and 6 is 
not too close to 90°. 
Throughout this paper, we work in the proton rest frame. 



nance condition dSOl l for the beam particles has two contribu- 
tions, namely the beam velocity and the offset from the beam 
center In some of the figures, we will plot "resonance lines", 
which represent the condition uji--kj^^\ = nfti. A solution of 
the resonance condition (Equation ( l30l l) corresponds to an in- 
tersection of a resonance line and a plot of a dispersion rela- 
tion Ui- = ujh- Since the distribution function of the alpha par- 
ticles has a finite thermal width, we also plot a spread of the 
resonance lines taking particles into account that are slightly 
faster or slower than the bulk speed of the beam. The two 
resonance lines with the same n then only represent the ap- 
proximate upper and lower limits of parallel particle speeds 
that are available in the distribution function. Generally, there 
are always particles that can fulfill the resonance condition in 
the area between both plotted resonance lines with equal n. 
This is the only point in our considerations where thermal ef- 
fects play a role, since we employ the cold-plasma dispersion 
relation. Some potential resonances are shown in Figure 2] 
Here we plot the six solutions in this wavenumber/frequency 
range to the dispersion relation. Equation (|3). The red lines 
show the first two (n = ±1) cyclotron resonances for alpha 
particles with parallel velocities V|| =U\± 0.05va. We take 



Ay I 



'II 



-U{\ to correspond to twice the parallel thermal 



speed Vthi = ^Ik^TwJmi of the alpha particles on the assump- 
tion that there are too few alpha particles with Ay|| > 2vthi to 
drive a significant instability. This situation corresponds to a 




c) 



p 




Vph 

Figure 5. Direction of the diffusive particle flux for cyclotron-resonant al- 
pha particles. The regions occupied by thermal protons (alpha particles) are 
shown as blue (red) filled circles. The aiTows and circles show the diffusion 
paths around the position of the parallel phase speed Vph for a cyclotron res- 
onance. The green dashed lines are iso-contours of the kinetic energy in the 
proton frame, a) U\ < I'ph, b) U[ > fph, c) Vpj, < 0. 

plasma with 



87rniA;Br||i / mnii 



Vthi 



10" 



(34) 



This means that the parallel particle velocity V|| in the reso- j^J^^ in Fisureffl 



where T^^ is the parallel ion (alpha-particle) temperature. The 
Landau resonance, which is only acting if Ei^, ^ or fc^ ^ 0, is 
located on resonance lines with a small spread around the line 
= f/;, which corresponds to the vicinity of the blue-dotted 



4. INSTABILITY THRESHOLDS 

In this paper, we focus on instabilities driven by the 
velocity-space diffusion of alpha particles that interact reso- 
nantly with the unstable wave. Some possible velocity-space 
diffusion paths for a beam of alpha particles are shown in Fig- 
ure |5] for the cyclotron-resonant case and in Figure |6]for the 
Landau-resonant case based on the assumption that the par- 
ticle distribution functions are isotropic about the mean ve- 
locity of the particle species. The resonant alpha parti- 
cles diffuse down the density gradient in velocity space, from 
larger particle concentrations towards smaller particle con- 
centrations along the diffusion paths. We concentrate on the 
normal case with positive wave energy first. In this case, if 
the particles gain energy during the resonant quasilinear dif- 
fusion, the energy is taken from the wave, and the wave is 
damped. On the other hand, if the particles diffuse towards 
smaller particle energies, then the energy lost by the particles 
is given to the wave, leading to wave amplification. To see 
whether particles gain or lose energy, the quasilinear-dif fusion 
paths can be compared with curves of constant particle energy, 
which are circles centered on the origin in the V||-v'j^ plane. 
The key point demonstrated by Figure |5] is that cyclotron- 
resonant wave-particle interactions cause alpha particles to 
lose energy if and only if < Vph < Ui. Therefore, if W > 0, 
the wavevectors and frequencies of unstable waves must lie 
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Figure 6. Direction of the diffusive particle flux for Landau-resonant al- 
pha particles. The regions occupied by thennal protons (alpha particles) are 
shown as blue (red) filled circles. The arrows show the diffusion paths at the 
position of the parallel phase speed Vpi, for a Landau resonance. The green 
dashed lines are iso-contours of the kinetic energy in the proton frame, a) 
J/, < Vph, b) J/, > Vph. 

above the line = and below the line Wr = k\\lJ\ in the k\y- 
plane. Figure |6] demonstrates that the same requirement 
is also valid for Landau-resonant particles. A mathematical 
proof of the instability condition < Vph < IJ\ is given in the 
Appendix. 

For simplicity, we take the protons to be Maxwellian for the 
purposes of inferring the effects of wave-particle interactions 
on the waves. We thus neglect proton temperature anisotropy 
and proton beams. As a consequence, if protons diffuse in 
velocity space from regions of large particle concentration to- 
wards regions of small particle concentration, then they must 
gain energy. Because the protons are the majority ion species, 
we expect that proton damping will dominate over the alpha- 
particle instability drive if thermal protons can resonate with 
a wave. 

The foregoing discussion and the condition cr, ^ 1 imply 
that a solution of the dispersion relation with k = ki, w^^i = uj\, 
and W > is driven unstable by resonant interactions with an 
alpha-particle beam that is isotropic about the beam velocity 
t/i in v-space if and only if the following conditions are satis- 
fied: 

1. The coordinates (k\^, wi) lie between two resonance 
lines with equal n with n = +1, n = 0, or n = -1. If 

uji) lies between two resonance lines with n = -1, then 
we say that the instability is driven by an n = -1 reso- 
nance, and likewise for « = and « = +1 . 

2. For an instability driven by an n = +1 resonance, E^^i 
must be nonzero (i.e., V ^+1). For an instability driven 
by an n = -1 resonance, Ek,i- must be nonzero (i.e., V ^ 
-1). For an instability driven by an « = resonance, Ej^^^ 
or must be nonzero. 

3. The parallel phase velocity of the waves is between zero 
and the beam speed: Q <uj\lk\\\ < Ui. 

4. The wave at ki and uji does not satisfy the resonance 
condition Equation dSOl l with thermal protons. 

Condition 2 above on the wave polarization follows from sub- 
stituting Equation ( |33] | into Equation ( |29] l. Importantly, a 
mode with W < that satisfies the above criteria is damped, 
because an increase in the energy of a negative-energy wave 
corresponds to a decrease in the wave's amplitude. In the next 
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Figure 7. Dispersion and polarization for U; = O.Si'a, t) = 0.05, and 9 = Q° . 
The color coding shows the polarization V. Only branch is shown. The 
two red lines are the resonance lines for n = —i. The dotted blue line shows 
the maximum frequency tj, = Ui, and the triple-dotted gray line is the min- 
imum frequency lU; = 0. 

two subsections, we use the above criteria to determine which 
branches of the dispersion relation are unstable and at which 
wavenumbers as a function of the beam velocity, the angle of 
propagation, and the parallel thermal widths of the alpha and 
proton distributions. 

4.1. Parallel Case (k± =0) 

The Landau resonance n = can be excluded for fast modes 
and cyclotron modes when k± =0 because E^, = (so that 
Landau damping is absent) and because transit-time damping 
(the other type of n = resonance) requires k± ^ ( Stix 1992). 
The counter-propagating alpha-cyclotron mode (branch C~ in 
Figure [U appears to be a candidate for an instability driven 
by the « = -1 resonance with the alpha particle beam, because 
it satisfies the four criteria for instability listed above when 
ojtr > 0, even for modest beam speeds such as Ui = 0.5va, 
which can be seen with the aid of Figure |7] However, this 
mode has negative energy when uj^r > 0, as illustrated in Fig- 
ure[3] and is therefore stable. 

If /3||i is sufficiently large, then alpha particles can sat- 
isfy the n = +1 resonance condition with the Alfven/proton- 
cyclotron wave (branch A in Figure [T) and potentially drive 
this mode unstable. To the best of our knowledge, such a 
parallel Alfven instability has not been described previously 
in the literature, although it might play an important role in 
the solar wind. We note, however, that our approximations 
break down when (3 > I, because we have approximated the 
wave frequencies using the cold-plasma dispersion relation. 
In order to investigate parallel Alfven/cyclotron instabilities 
at large rigorously, a more sophisticated approach based 
on the full hot-plasma dispersion relation would be needed. 

Although the counter-propagating alpha-cyclotron branch 
C~ is stable when k± = 0, the fast-magnetosonic/whistler 
mode becomes unstable at sufficiently large Ui due to the 
n = -1 resonance with the alpha particles. When Av|| = 
2vthi = 0.35va, the minimum value of Ui for this instability is 
~ 1.7va. With the given alpha-particle density, this situation 
corresponds to a plasma with /3||i ~ 0.006. These parameters 
are illustrated on the left-hand side in Figure |8] and it can be 
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seen from this figure that this instability satisfies all four crite- 
ria in the above list. This instability corres ponds to the paral- 
lel magneto sonic instab i lity dis cussed by i Daughton & Garvl 
(0998) and iGary et all (l2000al) . and sets in at a typical 
wavenumber of order the inverse proton inertial length Jlp/ v'a- 
The left-hand side in Figure [8] furthermore shows that at this 
drift speed but at a lower thermal speed the resonance lines 
do not intersect with the plot of the dispersion relation of 
the magnetosonic branch. This lower thermal speed repre- 
sents a situation with a lower provided the fractional num- 
ber density 77 is constant. As Av|| is decreased below 0.3va, 
it becomes increasingly difficult to excite this instability for 
two reasons. First, the value of U; needs to be increased 
simply in order to get the « = -1 resonance line to intersect 
the dispersion-relation curve for the forward-propagating fast- 
magnetosonic/whistler mode. Second, at these larger values 
of Ui, the fast-magnetosonic/whistler branch B merges with 
the counter-propagating alpha-cyclotron branch C~ through- 
out a finite interval of k^^ (i.e., they have the same value of Wii). 
Within this interval, W < for both modes, and any intersec- 
tion of an n = -1 resonance line with either mode's dispersion 
relation is unable to generate an instability according to the 
above four criteria. 

As Vthi is increased, our quasilinear/cold-plasma analy- 
sis suggests that the threshold for this instability moves to 
smaller Ui, because larger alpha-particle thermal speeds en- 
able the alpha-particle n = -1 resonance lines to intersect with 
the fast-magnetosonic/whistler dispersion relation at smaller 
values of Ui. This situation is shown on the right-hand side 
of Figure [8] With ?/ = 0.05, this situation corresponds to 
~ 0.02. The growth rate dependence of the magnetosonic 
instabili ty on the thermal speed o f the beam has been calcu- 
lated bv IDaughton & Garvl(ll998h for a proton-beam plasma. 
In qualitative agreement with our arguments, they have shown 
that the growth rate increases for higher temperatures of the 
beam. 

4.2. Oblique Case (k± =/0) 

A moderate obliquity does not significantly change the 
structure of the dispersion branches A, B, and in Fig- 
ure[T] It does, however, change the polarization. When k± y^O, 
the waves are in general elliptically polarized and can inter- 
act with particles via both the n = +1 and « = -1 cyclotron 
resonances. Furthermore, Landau-resonant interactions (with 
n = 0) are now possible, because Ei^^ becomes nonzero for 
oblique waves, and the nonzero value of k± allows for transit- 
time damping. 

The magnetoson ic in stability that was present in the par- 
allel case (Section 14.11 ) is also present in the oblique case. 
However, whereas the parallel magnetosonic mode resonantly 
interacts with the alpha particles only through the n = -1 cy- 
clotron resonance, the oblique magnetosonic mode can inter- 
act with alpha particles through either the n = -1 cyclotron 
resonance or the « = Landau resonance. For low drift speeds 
Ui > V'A, the resonance condition with n = is fulfilled at lower 
values for than the resonance condition with n = -l. The 
two resonance conditions are, therefore, expected to act on the 
magnetosonic mode in different wavenumber regimes. In the 
small-wavenumber limit, the fast-wave dispersion relation is 
oJkr = kv'A, and the instability condition U > uJkr/k\\ becomes 
Ui > va/cos0. This argument suggests that the threshold 
value of the drift speed for the obUque magnetosonic insta- 
bility is higher at larger 9. 



The intersection of the « = -1 resonance line with the dis- 
persion relation for proton-cyclotron waves can fulfill con- 
straint 2 for oblique waves due to their elliptical polarization. 
To fulfill constraint 4, the drift speed should be larger than 
Ui sa 1.2va. Then the intersection of the resonance line and 
the dispersion relation for oblique Alfven/proton-cyclotron 
waves occurs at low enough frequencies and wavenumbers 
(A:||Va/Op < 1) that proton damping can be ignored. This 
is the situation in Figure |9] which also shows the wave po- 
larization. We identify this instability with the Alfven I in- 
stability discussed by Gary et al. (2000b). It begins at U > 
L2va. There is a finite interval of k\\ in which the proton 
cyclotron wave and the counter-propagating alpha-cyclotron 
wave (branches A and C~ in Figure [T]i merge — i.e., in 
which the two waves have the same value of ujh-. Within this 
wavenumber interval, both waves have negative energy, and 
neither is driven unstable by resonant particlesfl 

Another instability that becomes possible at oblique prop- 
agation is the so-called Alfven II instability. We identify 
this instability with the forward-propagating proton-cyclotron 
mode, when this mode is driven unstable by Landau-resonant 
interactions with the alpha-particle beam. According to the 
third of the four criteria laid out at the beginning of Sec- 
tion ID this instability requires that Ui > LDp^^/k^. The parallel 
phase velocity LLiki/k\\ of the proton-cyclotron wave is va for 
^11 <C fip/vA and decreases towards zero at k\^\ 00. Thus, 
the condition Ui > ujkr/ ^[| is always satisfied for positive U 
provided that is sufficiently large. However, if is too 
large, criterion 4 in the above list is not satisfied, because 
ojkr ^ rip and the wave is damped by cyclotron resonant in- 
teractions with protons. In order to avoid proton cyclotron 
damping, must be < ilp/v'A (the precise upper limit de- 
pends upon the value of /3|[p for the protons). In order to have 
Ui > ijJkr/k\^\ at < fip/vA, U must exceed ~ 0.8va, which is 
the approximate threshold-value of U for this instability. This 
case is shown in Figure[TO] Criterion three in the above list of 
four criteria, combined with the Landau resonance condition 
^kr/k\\ = Vj|, imply that only alpha particles with V|| < U can 
drive the Alfven II instability. As U is increased beyond va, 
fewer and fewer alpha particles can undergo a Landau reso- 
nance with the proton-cyclotron wave, and thus there is an 
upper limit on U, for the Alfven II instability. As increases 
for the alpha particles, this upper limit increases. 

There is evidence in previous treatments of beam insta- 
bilities that another kinetic effect may weaken or suppress 
the oblique instabili ties in certain para meter regimes. It 
was found by Mont gomery et al.l (119761) in a proton-beam 
plasma that the oblique Alfvenic instability (which corre- 
sponds to Alfven II according to Daughto n & Garvl IT998) 
shows a higher growth rate at greater values of Ti/Tp, whereas 
the cyclotron-resonant m agnetosonic instability does not. 
iMontgomery et all d 1 9761) interpreted this result as a conse- 
quence of Landau-damping of these modes by electrons. At 
higher electron temperatures, the gradients of the electron 
distribution function in velocity space become smaller since 
the width of the distribution function increases. The influ- 
ence of electron Landau damping is determined by these gra- 
dients and, therefore, mitigated at higher electron temper- 
atures. In our model, the Landau-resonance condition for 

Within this wavenumber interval, the wave frequency is complex, and for 
one of the two branches the solution to the cold-plasma dispersion relation has 
a positive imaginary part. We neglect this type of instability in our discussion 
for the reasons discussed in the Introduction. 
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Figure 8. Dispersion relation and polaiization of the branches B and C for t\ = 0.05 and 9 = 0. The color coding shows the polarization V. The resonance lines 
for n = —1 are shown for different thermal speeds. The dotted blue line shows the maximum frequency tOr = k^^ U{, and the triple-dotted gray line is the minimum 
frequency cjr = 0. Left: U\ = I.Vva- The thermal spread corresponds to Av|| = 2vii,i = O.OSi'a for the red resonance lines and to Av|| = 2vthi = 0.35va for the green 
resonance lines. Right: U\ = 1.5va and Ai'ii = 2vthi = 0.7va. 
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Figure 9. Dispersion relation and polarization for the wave branches A and 
C". The same parameter set as in Figure|4]is used. The color coding shows 
the value of the polarization function V. 

electrons would correspond to resonance lines that are very 
widely spread due to their large thermal speed. Therefore, 
almost all waves with ^ in our frequency range can un- 
dergo electron Landau damping. The Alfven I instability may 
not be as strongly influenced by this effect as the Landau- 
resonant instabilities since the Alfven I instability generally 
shows higher growth rates in the numerical solutions of the 
hot-plasma dispersion relation under typical parameters (e.g., 
Garv et al.l' 2000b.). This effect may also be the reason why 
Montgome ry et al.l (119761) find that the oblique magnetosonic 
instability has a significant growth rate over only a very small 
range of the plasma parameters, and why Gary et al. (2000b) 
do not discuss this mode for a proton-alpha plasma. The de- 
tails of this effect, however, require a deeper study and are 
beyond the scope of this work. 

5. DISCUSSION AND CONCLUSIONS 



Figure 10. Dispersion relation and polarization for the wave branch A with 
8 = 45°, rj = 0.05, and U\ = O.Si'a- The color coding shows the value of the 
polarization function P. The resonance lines represent the Landau resonance 
with n = and Avy = 2V(hi = O.Ii'a. Only resonance lines below the line 
u)s = k\\U\ can drive the Alfven II instability according to criterion 3. 

In this paper, we use quasilinear theory to explain 
the thresholds for Alfven/ion-cyclotron and magne- 
tosonic/whistler instabilities driven by resonant interactions 
with an alpha-particle beam in low-/? plasmas. Our arguments 
rely upon three properties of the wave modes in question: 
the real part of the wave frequency Wjtr, the polarization of 
the wave electric fiel d, an d the sign of the wave energy den- 
sity W . (See Section |23] for a discussion of negative-energy 
waves.) To obtain approximate values for these quantities as 
functions of the wavevector k, we have used the cold-plasma 
approximation, taking into account the finite electron mass. 

The general considerations that determine whether a wave 
is driven unstable by resonant interactions with a beam of al- 
pha particles are the following. First, in order for a wave to 
be unstable at some value of kn, it can not undergo cyclotron 
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resonant interactions with the thermal protons at that k\\ , since 
such interactions would strongly damp the wave. Second, the 
wave and alpha particles must satisfy the wave-particle res- 
onance condition. Equation ( l30l ). Third, the resonant wave- 
particle interactions must cause the alpha particles to lose en- 
ergy, so that the wave can gain energy and grow. (This state- 
ment assumes that W >0 — more details on the stability and 
instability of negative-energy waves are given in Section |4]l. 
As we show in Section|4]and in the Appendix, resonant wave- 
particle interactions cause alpha particles with a distribution 
function that is isotropic about their mean velocity to lose en- 
ergy when < LOui/kw < Ui, where Ui is the speed of the alpha 
particle beam, and all quantities (e.g., w^,, U;) are measured in 
the proton rest frame. Fourth, the wave polarization must sat- 
isfy the following conditions. In order for a wave to be driven 
unstable by an n = +1 (n = -1) cyclotron resonance, /s^j (/s^.r) 
must be nonzero. When k± = 0,a wave can be driven unstable 
by the n = Landau resonance only when E^, ^ 0. 

When we apply these considerations to parallel waves with 
= 0, we find that alpha particles are unable to drive reso- 
nant instabilities when Av|| = 2vthi is much smaller than 0.35. 
When Av|| ~ 0.35, our criteria suggest that the magne- 
tosonic/whistler wave becomes unstable via the n = -1 res- 
onance when Ui exceeds ^ 1.7va, consistent with pre vious 
numerical results (iLi & Habball2000HGarv et al.ll2000al ). Our 
arguments suggest that this minimum threshold on Ui de- 
creases slightly as the parallel thermal speed of the ions in- 
creases further Our arguments also suggest that as Vthi (and 
therefore /3||i) increases, the alpha particles eventually drive 
the proton-cyclotron wave (branch A in Figure [T) unstable in 
the n = +1 resonance. We note, however, that our conclusions 
about the high-/3||i case must be treated with caution, because 
of our use of the cold-plasma dispersion relation. 

When we apply the above instability criteria to oblique 
waves with k± ^ 0, we find that the magnetosonic mode 
can still be unstable, but that it can now be rooted in ei- 
ther the cyclotron resonance (n = -1) or the Landau reso- 
nance (n = 0). The minimum threshold to excite the oblique 
magnetosonic instability via the Landau resonance is some- 
what smaller than for the parallel magnetosonic instability, 
and also less dependent on VtM. The wavenumbers at which 
the mode is unstable when driven by the n = resonance are 
smaller than the unstable wavenumbers for oblique magne- 
tosonic instabilities driven by the n = -1 resonance. At these 
smaller wavenumbers, the dispersion relation of the magne- 
tosonic mode for the n = resonant instability can be ap- 
proximated by ojkr ~ k^\VA/ cos 9 which leads to an instability 
threshold of U ^ v\/ cosO according to the instability crite- 
rion 3 introduced in Section |4] At oblique propagation, there 
are also two types of Alfven/cyclotron instabilities that are 
not present at parallel propagation. The first of these, the 
"Alfven I" instability, is driven by the « = -1 cyclotron res- 
onance. Our analysis suggests that this instability requires 
U > 1-2va , which is consistent with the numerical results of 
iGarv et alT(l2000bl) for a plasma with Tg/Tp = 4. Our physical 
interpretation of this instability threshold is that it is the min- 
imum value of U for which alpha particles can resonate with 
the proton-cyclotron wave (branch A in Figure [T} through 
the n = -1 resonance at a wavenumber that is sufficiently 
small to avoid strong cyclotron damping by the protons. The 
second obhque i nstabilit y, the so-called Alf ven 11 instability 
dPaughton & Ga rv 1998; Gary et al.ll2000bl) . is driven by the 
n = Landau resonance. Our arguments suggest that this 



Table 1 

Summary of the Unstable Modes and Their Thresholds 



Name 


Resonance 


Threshold 


Parallel magnetosonic 


«=-l 


~ 1.7VA 


Oblique magnetosonic" 


« = or H = -1 


> I'A / cos 6 


Alfven I 


n = -l 


~ 1.2VA 


Alfven II 


n = 


~ 0.8vA 



Note. — The three obliqu e ins tabilities may depend on 
Te/7p as described in Section|4]2] 



"The oblique magnetosonic mode can be driven unstable 
by two different resonances. The n = resonance has the 
lower threshold given in this table and occurs then at lower 
wavenumbers compared to the n = —l resonance. 

instability require s that U > 0.8va, consistent with the nu- 
merical results of IGarv et al.l (I2000bl) for the case in which 
Te/Tp = 4. This lower threshold on U is the minimum beam 
speed for which the proton-cyclotron wave (branch A in Fig- 
ure [Til can undergo a Landau (« = 0) resonance with the al- 
pha particles at a that is sufficiently small to avoid strong 
proton cyclotron damping. We note, however, as discussed 
in Section l4~2l that the Alfvenic instabilities and the oblique 
magnetosonic instability may weaken as Te/Tp decreases be- 
cause of the more efficient electron Landau damping of these 
modes. This effect may lead to a weakening or suppression 
of these instabilities in the fast solar wind, where is typi- 
cally smaller than Tp by a factor of a few. A summary of our 
findings for the instability thresholds of oblique and parallel 
waves IS given m Table m 

The principal utility of our analysis is that it provides a con- 
ceptual framework for understanding previous numerical re- 
sults on the different instability thresholds, including the find- 
ing that the instability threshold drops from U > 1.7va for 
parallel propagation to Ui > O.Sv'a for oblique propagation. 
Our discussion also describes the different conditions under 
which cyclotron resonances (either n = +l or« = -l) and the 
Landau resonance n = are important. 

The wider importance of this topic is that beam-driven 
instabilities place an interplanetary 'speed limit' on alpha- 
particle beams. As demon strated in numero us previous stud- 
ies (Marsch & Li vi 1987;: Ipaughton et al.1 [1999: Gary et all 
2000b; Golds teinetaiT 120001; iReisenfeld et all 120011; 
Kaghashvili et al.' 120031; 'Gary et alj 120031; iHelUnger et alj 
2003; Hellinger & Travm'cek 201 1|K when U exceeds the 
threshold for a beam instability, exponential growth of the 
unstable wave and the resulting wave-particle interactions 
decelerate the beam and limit the beam speed to approxi- 
mately the marginally stable value. If this scenario is correct, 
then it is the instability with the smallest t/i-threshold that is 
most important for limiting the beam speed in the solar wind. 
At low /3, the instability with the lowest threshold is the 
Alfven 11 instability, at le ast w hen T^/Tp is sufficiently large 
(see discussion in Section |4~2] |. At Tf^/Tp < 1, the Alfven II 
instability may be more difficult to excite, in which case the 
instability with the lowest instability threshold could become 
the Alfven I instability. Regardless of which instability has 
the lowest beam-speed threshold, each of the instabilities 
that we have investigated has a threshold-value of U that is 
normalized to the Alfven speed. Because va decreases with 
increasing heliocentric distance r outside of the solar corona, 
if instabilities limit the value of U to ~ va in the solar wind, 
then they lead to ongoing deceleration of the alpha-particle 
beam at r > rmin, where r^in is the radius at which U first 
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reaches a value ^ va as the alpha particles flow away from 
the Sun. This scenario of deceleration by beam-driven 
instabilities is consistent with Helios measurements showing 
that Ui < Va at 0.3 AU < r < 1 AU and Ulysses observations 
of proton beams at up to r = 3 AU jMarsch et al.l 119821; 
iGoldstein et al. 2000). 

Another reason that beam-driven instabilities could be im- 
portant in the solar wind is that they could lead to an inverse 
cascade of wave power from large k\\ to smaller k\\. In this 
picture, the ion beams generate Alfven/cyclotron waves at 
cj ^ ilp and ^ rip/v'A. In the MHD Hmit {k ^ VL^/va), 
an Alfven wave can decay into another Alfven wave, which 
propagates in the opposite direction, and a slow mode wave, 
which propagates in the same direction as the mother wave, a 
process called "the parametric instability" (ISagdeev & GaleevI 
[1969; Cohen & Dewar 1974). The daughter Alfven wave has 
a slightly smaller k\\ and u) than the mother wave, and over 
time the par ametric instability leads t o an inverse cascad e, as 
described bvlCohen & Dewarl (119741) . iKuznetsovl (1200 Ih . and 
iChandranI dSoOSl) " for a weak turbulence scenario at low beta 
and k <C rip/vA. Hybrid numerical simulations have shown 
that a similar parametric instability and inverse cascade occur 
at k ~ r^n/vA (M arkovskii et al. 2009; Matteini et al. 2010; 
IVerscharen et al.ll2012bl) . An inverse cascade at A; ^ i^p/vA 
feeds energy into the low-frequency MHD regime, where 
the MHD inverse cascade can continue transferring energy 
to lower frequencies. This scen ario is c onsistent with obser- 
vations from the ACE satellite (IStawar z et al. 2010), which 
show evidence for energy transfer from small to large scales 
when the solar wind is in a state of high cross-helicity. An 
inverse cascade could produce "slab-like" magnetic fluctu- 
ations in the solar wind with k\\ < k± over a broad range 
of parallel wavenumbers, w hich could lead to efficient scat- 
tering of energetic particles (Jokipii" 1966; K ulsrud & Pearcd 
[1969). In contrast, the quasi two-dimensional fluctuations 
with k± ^ k \\ produced by the f orward Alfven- wave energy 
cascade (iShebaU n et al.' 498 1 iGoldreich & SridhaJ [19951 
lOuphton etalT Tr995; Verschar en et al. 1 120 12ab are high ly m- 
effective at scattering energetic particles (iBieber et al.lll994t 
[Chandran 2000; Yan & Lazarian 2002). Further investiga- 
tions of this inverse cascade will be an interesting task for 
a future study, which could eventually connect plasma insta- 
bilities and turbulence to the transport of energetic particles in 
the solar wind. 
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APPENDIX 

proof that an ISOTROPIC DISTRIBUTION OF BEAM 
PARTICLES LOSES ENERGY FROM RESONANT 
INTERACTIONS WITH WAVES SATISFYING < Wfa/zty < f/, 

In this Appendix, we consider the way that the particle ki- 
netic energy changes in response to resonant wave-particle in- 
teractions. The kinetic energy density in the proton rest frame 
of particle species s is £ = 0.5ms J fs\^d^v. Since we have 
taken /, to be independent of spatial position x, £ is a func- 
tion of t alone, and we can write 



d£ 1 
df 2 



2 ^/i j3 



Using Equation (l27l i to describe the quasilinear evolution 
of /;, we can rewrite Equation (lAlt as 



df 
dT 



1™ d^/t / dvii 



X / dv_L 



+°° 2,2 



V 



Gv 



xd(ujkr-k\iv\i-nns)\'>pn.k\ Gfs. (A2) 

We assume that /, is isotropic about the beam velocity Uitz 
and monotonically decreasing with distance from UiC,^ in ve- 
locity space, i.e. 



/, = /,((v-f/,e,)2), 



(A3) 
(A4) 



where the prime indicates t he d erivative with respect to the 
argument of fs in Equation (IA3l l. Applying the G-operator to 
fs yields 



G/, = 2yi( 1-^ )/;. 



(A5) 



The first G-operator in Equation (IA2t can be eliminated by 
integration by parts, leading to 



df , 



(A6) 



where 



4=1 



k\\Ui 



1. 



1 +00 



lim 



X / dvi^(5(a;fa--^||V||-«rjj)|V'n,iPl/i| 



(A7) 



The quantity 4 describes the contribution to d£/df from 

waves at wave vector k. 

All terms on the right-hand side of Equation (IA7) are posi- 
tive semi-definite except for the quantity 



'= 1- 



(AS) 



which can be positive or negative depending on the ratio of 
U\ to Vph = W(tr/^|| and the sign of Vph. The sign of 4 is then 
negative if and only if: (1) g < 0; and (2) |V'«.jtP 7^0 and 
f[ y (meaning /' < given Equation ( IA4I) ) for at least some 
subset of the particles that satisfy the resonance condition. 
The case in which f[ = for all resonant particles is not 
relevant to ion beams in the solar wind, and so we do not 
consider this case further The quantity iV'n.jtp can vanish 
if the wave amplitude vanishes or if the wave polarization 
and propagation direction satisfy certain constraints: e.g., if 
« = 1, A:^ = 0, and E/^ i = 0. We have addressed the constraints 
on wave polarization elsewhere in this paper (see, e.g., 
instability criterion 2 in Section [4), and we assume that 
there are waves that can resonate with the particles. In the 
remainder of this Appendix, we thus neglect the case in 
(Al) which IV'n.itP is identically zero in the integrand within the 
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brackets of Equation ( IA7b . The sign of 4 is then equal to 
the sign of g. Therefore, if we wish to determine whether 
waves at some particular wave vector ki act to increase or 
decrease the particle kinetic energy, we need only determine 
the sign of § at k = kj , which depends only upon the value 
of ujty/ at k = ki . The condition ^ < is equivalent to our 
criterion 3 in Section |4] namely < utitr/k^^ < Ui. We note 
that if a spectrum of waves is excited over a range of k values, 
then some waves may act to increase £ while others act to 
decrease £. The sign of d£/dt m ust th en be determined by 
evaluating the integral in Equation (IA6b . 
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